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Abstract 

Kar's recent proof showing that a maximally entangled state of 
two spin-1/2 particles gives the largest violation of a Bell inequality is 
extended to N spin-1/2 particles (N > 3). In particular, it is shown 
that all the states yielding a direct contradiction with the assumption 
of local realism do generally consist of a superposition of maximally 
entangled states. 



Recently, Kar (see Ref. Q, and references therein) has shown that a 
maximally entangled states of two spin-1/2 particles not only gives a max- 
imum violation of the CHSH inequality ||] but also gives the largest vio- 
lation attainable for any pairs of four spin observables, these pairs being 
noncommuting for both systems. To prove this, Kar made use of an elegant 
(and powerful) technique based on the determination of the eigenvectors and 
eigenvalues of the associated Bell operator Q. In this Letter we would like 
to extend these results to the case in which N spin-1/2 particles (N > 3) are 
considered. We will show that the most general iV-particle state giving the 
largest violation of a Bell inequality does consist of a superposition of maxi- 
mally entangled states. As expected, those states giving maximal departure 
from classical expectations correspond to the class of states introduced by 
Greenberger, Home, and Zeilinger in proving Bell's theorem without using 
inequalities |f]. [5| . 

In order to look for a violation of local realism when dealing with N 
spin-1/2 particles it is necessary to consider correlation functions involving 

*An almost identical version of this paper was originally published in: J.L. Cereceda, 
Phys. Lett. A 212, 123 (1996). 
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measurements on each of the particles. A suitable generalisation of the 
CHSH inequality to an Admeasurement scheme was obtained by Hardy ||. 
Hardy's inequalities can be written in the form 

- 2 < (B n ) < 2 , (1) 

where (-E>h) denotes the expectation value of some Hermitian operator I?h 
(the so-called Bell operator [||) acting on the 2 Ar -dimensional tensor prod- 
uct space associated with N spin-1/2 particles. For concreteness, and for 
purposes of comparison with Ref . Q , in what follows we concentrate on the 
case N = 3. Later on, we shall consider the case of arbitrary N. For three 
spin-1/2 particles the representative Bell operator is || 

B H = a (ni) a (n 2 ) a (n' 3 ) + a (n[) a (n 2 ) a (n' 3 ) 

+ ° ( n 'i) a ( n 2) <* ( n 3) - cr (ni) a (n 2 ) a (n 3 ) , (2) 

where the observable operator a (n^) (a (n-)) corresponds to a spin measure- 
ment on particle i (i = 1,2,3) along the direction nj (n-). The square of 
this Bell operator is given by 

Bl = 41- [aW^K)] [a(n 2 ),a(n 2 )] 

- [a- (n 2 ) , a (n' 2 )] [a (n 3 ) , a (n 3 )] 

- [a (m) , a (ni)] [a (n 3 ) , <x (n' 3 )] . (3) 

Now, replacing each commutator [a (nj) , a (n-)] by its valueQ 2i sin0j <7j_j, 
we get 

Bft = 41+4 sin 6<i sin0 2 c±i °J_2 
+ 4 sin 2 sin 3 crj_ 2 crj_ 3 

+ 4 sin 0i sin 3 a±i crj_ 3 , (4) 

where 0j is the angle included between nj and n-, and a±i is the spin operator 
for particle i along the direction perpendicular to both nj and n-. From 
expression (4) we can see at once that the largest eigenvalue of B^ is 

fMi = 4 (1 + I sin 0i sin0 2 | + |sin0 2 sin0 3 | + |sin0i sin0 3 | ) , (5) 

which attains a maximum value of /i max = 16. From (4) it is also apparent 
that to every eigenvalue of B^ there corresponds a pair of (degenerate) eigen- 
vectors, namely, |ct_li, ctj_ 2 , a±s) and | — a±i, —a±2, —^±3), where \ a±i) 

1 This expression for the commutator actually differs from that of Ref. jl]] by a factor 
— 1. This is because, following the usual convention, we take here the angle 9i as the 
negative of that used in Ref. Jl[. 
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(I — a ±i)) is the eigenvector of a±i with eigenvalue a±i = ±1 (— a±i = =Fl). 
In particular, the eigenvectors corresponding to the largest eigenvalue (5) are 
(in an obvious notation): | j, f, | ) and ||, j, J.) for sgn (sin#i) = sgn (sin 62) 
= sgn(sin6> 3 ); ||, T, I) and ||, |, T) for sgn(sin6»i) = sgn(sin6> 2 ) 
/ sgn(sin6» 3 ); ||, I, I) and ||, |, f) for sgn(sin6» 1 ) ^ sgn(sin6> 2 ) 
= sgn(sin6'3); and, finally, |T, I, T) and ||, |, [) for sgn(sin6q) 
= sgn (sin #3) 7^ sgn (sin # 2 ). On the other hand, it can be easily seen that 
the minimum possible eigenvalue for B# is zero. So, for example, |f, T> I) 
is an eigenvector of B^ with zero eigenvalue whenever 6\ = 62 = #3 = vr/2. 
Of course, the operator B^ cannot have negative eigenvalues because this 
would imply the (Hermitian) operator Bn has a complex spectrum. 

As every eigenvalue for 5 H must lie in the interval [0, 16] it follows that 
the eigenvalues for I?h are necessarily restricted to lie in the interval [—4, 4]. 
Consequently, inequality (1) for the Bell operator (2) will be violated for 
those eigenvectors of -Bh with eigenvalues A fulfilling 2 < |A| < 4. Note 
that the maximum amount of violation of Hardy's inequality predicted by 
quantum mechanics is by a factor of 2, instead of the factor \[2 achieved 
in the CHSH inequality. Moreover, it is worth pointing out that the results 
mentioned in the present paragraph remain valid for an arbitrary number of 
particles. This is so because Hardy-type inequalities involve only four cor- 
relation functions regardless of the number N of measurements considered 
IJ. This fact will be used later in considering the A^-particle case. 

In view of Eq. (5), the largest eigenvalue of -Bh will be 

|A/| = 2(1 + |sin0i sin(9 2 | + |sin(9 2 sin<9 3 | + |sin6q sin 6» 3 | ) 1/2 . (6) 

Now, as a product state cannot give rise to violations of local realism (see 
below), it follows that an eigenvector of Br with eigenvalue A such that 
2 < |A| < 4 must necessarily consist of an entanglement of the two degener- 
ate eigenvectors of -B H with eigenvalue A 2 ||| , 

\^) = a \zx, z 2 , z 3 ) + (3e^\ -z\, -z 2 , -z 3 ) , (7) 

where, for simplicity, we restrict ourselves to directions rij and lying in the 
x-y plane (these directions being specified by the azimuthal angles fa and 
respectively), so that \zi) ( | — Zi)) designates the eigenvector of the spin 
operator along the z-axis for particle i, with eigenvalue Zi = ±1 (— Z( = ^1). 
It turns out (see Eq. (8) below) that the relative signs of z%, Z2, and z% 
in Eq. (7) are uniquely determined by A (for fixed values of 9±, 62, and 
#3.) Likewise, the real coefficients a and (3 (which are assumed to satisfy 
the normalisation condition a 2 + /3 2 = 1), as well as the phase factor <p 
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will depend on the eigenvalue A. We will now show that for the case in 
which the eigenvector (7) is associated with the largest eigenvalue (6), these 
coefficients must fulfil the condition \a\ = \(3\ = l/y/2. In other words, the 
largest possible violation occurs for maximally entangled states. This can 
be seen by directly evaluating the expectation value for the state vector (7). 
This is given by 

(* \B H \ *) = 2 a/3 [ cos (0 - z^ - z 2 <p' 2 - ^3) 

+ COS - Z\4>\ - Z 2 fa ~ Z 3 fa 3 ) 
+ COS ((/) - Z X fa x ~ Z 2 fa 2 ~ z 3<fe) 
- COS ((f) - Z\fa ~ Z 2 (f>2 ~ Z 3 fa)] ■ (8) 

As the product a(3 factorises out in this expression, it is clear that in order for 
(8) to reach its largest value (6) it is necessary that the absolute value of a(5 
be a maximum, i.e., \a\ = \(5\ = Although it is apparent from (8) that 

this must be the case for |A/| = 4, it might seem to be worthwhile checking 
explicitly the above statement for the case that |A/| < 4 by considering a 
concrete example. So, let us take the values fa = fa = 03 = 0, fai = fa 2 = 
7r/2, and fa 3 = 7r/4, so that 9\ = #2 = Tr/2 and #3 = 7r/4. For these values we 
find that (see Eq. (6)) |A/| = 2 (2 + \/2 ) 1/2 . On the other hand, the absolute 
value of (*|B H |^) is found to be, 4\a(3\ | 2" 1 / 2 sin - (1 + 2~ l l 2 ) cos</> |, 
where we have put z\ = Z2 = £3 = +1 in Eq. (8) since, for the above 
values for 6i, we have sgn (sin#i) = sgn (sin #2) = sg n (sin #3) (of course the 
reasoning remains essentially unchanged if we instead choose z\ = Z2 = Z3 = 
— 1; the important point is that sgnzi = sgnz2 = sgn Z3). Therefore, for 
some a, (3, and (ft, we must have 

2\a(3\ 2- 1/2 sin0- (1 + 2~ 1/2 ) cos = (2 + v^) 1 ^ . (9) 

The only values for which this equality holds are \a\ = = I/a/2, and 
4> = —it/8 + n7r, where n = 0, ±1, ±2, ... . This follows at once from the 
fact that the function /i(</3) = | 2 -1 / 2 sin </3 — (1 + 2 -1 / 2 ) cos <f> \ reaches its 
maximum value for 4> = —ir/8 + nir, and this value is precisely (2 + \/2) 1//2 - 
From the preceding example it is obvious that, if expression (8) is to be 
equal to the largest eigenvalue (6), the phase factor <p must be a suitable 
function depending only on the angles 4>i an d <p[ (or, equivalently, on the 
angles fa and 0j, as 0j = fa { — fa). This is so because for the largest eigenvalue 
case the condition 2 \a(3\ = 1 is always met, and then Eq. (8) involves only 
the variables (ft, (pi, and fa. This dependence can be easily obtained for 
the special (and important) case where 9i = n/2 (i = 1,2,3). In this case 
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the eigenvalue (6) attains its maximum value 4, and then the following four 
equalities should be simultaneously fulfilled 



COS (( 


P - 4>\ - ■ 


02 


-0' 3 ) 


= ±1 


COS (( 


i> -</>[-■ 


4>2 




= ±1 


COS (( 


i> -</>[-■ 


4>2 


-03) 


= ±1 


COS (( 


i> - <t>i - ■ 


4>2 


-<h) 


= Tl 



(10a) 
(10b) 
(10c) 
(lOd) 

where we have put Z\ = Z2 = z% = +1 in (8) since, as before, sgn (sin6q) = 
sgn(sin#2) = sgn (sin 6*3). So, recalling the relationship = <pi + vr/2, it is 
a trivial matter to see that equalities (10a)-(10d) can be matched provided 
that <j) = (f)\ + (f>2 + 03 + mr, n = 0, ±1, ±2, ... . It thus follows that the 
state vector (see Eq. (7)) 

|¥> = -L ( I +,+,+> ± ^+<h+<h) 1 _ _)) f 

will be an eigenvector of the Bell operator 

B H = a a{<t>i) a(0 2 + vr/2) a(0 3 + vr/2) + (3 a(0i + vr/2) <r(<h) ^3 + tt/2) 
+ 7 <t(0i + tt/2) a{4> 2 + tt/2) <t(0 3 ) + 8 o(<t>l) o(<h) <r(<h) , 

(12) 

with eigenvalue =f4, whenever a = (3 = 7 = — <5 = +1. More generally, it 
can be shown that any state of the form 



I* 



) = 4= ( I zi, *2, z 3 ) ± ^1+^+^3) 1 _ Z2; _ Z3)) 5 



(13) 



is an eigenvector of the Bell operator (12) with eigenvalue +4 or —4 for a 
suitable choice of the sign factors a, [3, 7, and 5 (provided, in any case, 
that af3^5 = — 1). In this way, it is clear that the product of the quantum 
expectation values for the operators T% = a{4>\) a{(j>2 + 7r/2) cr(0 3 + 7r/2), 
T 2 = <r(fo + vr/2) afa) <r(fa + tt/2), T 3 = (r(0 a + tt/2) <t(0 2 + vr/2) <7(0 3 ), 
and T4 = cr(^>i) c(</>2) &{<j>z), must be equal to —1 when evaluated for any of 
the states (13). Indeed, the set of vectors (13) forms a basis set (with a total 
of eight linearly independent vectors) which simultaneously diagonalises the 
four (commuting) operators 7\, T 2 , T 3 , and T4 J7|]. As a result, the value 
— 1 for the above product of expectation values actually arises from the fact 
that Ti T2 T3 T4 = —I. That these quantum predictions for the expectation 
values indeed lead to a direct contradiction with the assumption of local 
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realism constitutes the theorem of Greenberger, Home, and Zeilinger [Q, 
(see also Refs. Q ||.) In fact, it can be easily shown j|] that a maximum 
violation of inequality (1) always entails a nonlocality contradiction of the 
GHZ type. 

Before analysing the iV-particle case, it is worth noting an implication 
of Eq. (8). Indeed, from that equation it follows that, for fixed (ft, the eigen- 
value A does not determine the individual values of a and (3 but, rather, 
the value of their product a(3. This means that, whenever \a\ ^ \/3\, if the 
state vector Y& a p) = a\z\, z%, z%) + (3e l ^ \ —z\, —Z2, —z%) happens to be 
an eigenvector of (2) with associated eigenvalue A, the same is true for the 
(linearly independent) vector \^f/3 a ) = P z%, z%) + a e 1 ^ \ —z\, —Z2, —z^). 
This degeneracy is due to the very structure of the state vector (7), and 
will be called here a trivial degeneracy. Notice that the trivial degeneracy 
is removed when A corresponds to the largest eigenvalue since, as we have 
seen, in this case we have \a\ = \j3\ = 1/V2. This type of degeneracy is to 
be distinguished from the nontrivial degeneracy which occurs when eigen- 
vectors with different relative signs for the z's are associated with the same 
eigenvalue. As already noted, for |A| > 2 the relative signs of z±, z%, and zj, 
in (7) are uniquely determined by A so that any eigenvector of the Bell op- 
erator (2) violating inequality (1) (excepting the nondegenerate eigenvector 
corresponding to the largest eigenvalue) is only trivially degenerate. 

Turning to the A r -particle case, one could equally prove that an N- 
particle state giving the largest violation of Bell's inequality has to be maxi- 
mally entangled. Properly speaking, such an A r -particle state will in general 
consist of a superposition of maximally entangled states. That this requeri- 
ment has to be met for those states yielding the maximum violation follows 
in a rather straightforward way when one considers the correlation function 
P{.4>ii ^2, • • • j <t>N\ ^f) that quantum mechanics predicts for a general (pure) 
state of the form 

1^) = a z u z 2 ,...,z N I 21, Z 2 ,... , Z N ) , (14) 

21, 22,... , 2JV 

with 

n —1(7 I J0z 1 ,z 2 ,...,z N ST^ I I 2 — 1 

" 21, 22,... , 2JV I U 21, 22,... , 2JV I C , / j | " 21, 22,... , 2]V I 1 , 

21, 22,... , Z N 

and where \zi) represents the eigenvector of the spin operator along the z-axis 
for the ith particle (i = 1, 2, . . . , N), with eigenvalue Zi = ±1. As before, 
and without loss of generality, we assume that each particle is subjected to 
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a spin measurement along a direction lying in the x—y plane, with azimuthal 
angle fa. The quantum prediction for P(fa, fa, ■ ■ ■ , <Pn] ^) is given by 



P(fa, fa, . . . , 4> N ; ^/) = 2 J2* \ a z 1> z2,...,z N \\a-g u - Z2 ,...,-z N \ 

21, Z 2 ,... , Z N 

x cos((j) Zl}Z2v .. >ZN - Zlfa - Z 2 fa Z N (f) N ) , 

(15) 

where 4> Zl ,z 2 ,...,z N = -z u -z 2 ,...,-z N ~ & zi,z 2 ,...,z N , and where the asterisk 
indicates that, for each pair of combinations of indices z\, z 2 , ■ ■ ■ , Zjv and 
—zi, —Z2,... , —zn, the summation runs over either the indices z\, z 2 , - ■ ■ , 
zn or —z%, —Z2, ■ ■ ■ , —zjq. By examining Eq. (15), one finds that for the 
function P(fa, fa, ■ ■ ■ , <Pn] to take on its extreme values ± it is necessary 
that \a Zl ,z 2 ,...,z N \ = I °-zi, -z 2> ... , -z N \ for all z\, z 2 , . . . , £jv- This is because, 
if this condition is fulfilled, then we have 2 Z2 z | a Zl ,z 2> ... )ZN \ 2 = 1> 
and thus expression (15) may attain the value ±1 for a suitable choice of the 

relative phases <p Zu z 2 ,...,z N , namely, for (j> zl ,z 2 ,...,z N = zifa + z 2 fa H h 

zn<Pn + n7r , where n = 0, ±1, ±2, ... . It thus follows that a state vector 
of the form 

l*>= C z u z 2 ,...,z N \^(zi,z 2 ,... , z N )) , (16) 

21, Z 2 ,... , Z N 

with Ht uZ2 ,...,z N \cz 1 ,z 2 ,...,z N \ 2 = 1, and 

\V(z h Z 2 ,..., Z N )) = -1= ( I Zl , Z 2 , . . . , Z N ) ± e i(-i^+^^+-+^^) 

x I -zi, -z 2 , • • • , -2jv) ) , (17) 

will yield the value ±1 for P(fa, fa, ■ ■ ■ , 4>n', *&), and then it will be able 
to violate maximally a Bell inequality built up from correlation functions 
involving spin measurements in the x-y plane. State vector (13) can be 
regarded as the simplest instance of Eq. (16). In this case only one term 
appears because either of the eigenvectors associated with the extreme eigen- 
values A = 4 or A = —4 of the Bell operator (12) is nondegenerate. In general 
the number of terms appearing in expansion (16) will be equal to the di- 
mensionality of the (nontrivially) degenerate subspace corresponding to the 
maximum eigenvalue of the relevant Bell operator. 

The Bell operator Br for the Hardy inequality (1) takes the general form 
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B n = a a (0?) a • • • a + f3 a (fi) a (4) ---a (<) 

+ 7 o (#) a ($) • ■ ■ a (<&) + 5 a (cf>{) a • • • a (<) , 

(18) 

where, as before, a, (5, 7, and 5 are sign factors with af3j5 = — 1. As was 
mentioned, the state vectors violating inequality (1) will be those eigenvec- 
tors of (18) with eigenvalues A such that 2 < |A| < 4. Naturally, the four 
parameters (f)f , </>f, <j)] , and 4 are not all independent. Indeed, for each 
value of i, there are the following possibilities ||: (i) 4 = <f>f , 4 = 4>]; 
(ii) (/)] = <fif, </)f = 4; and (iii) (pf = 4>f , <f>] = 41. ■ I n an Y case the most 
general eigenvector for the Bell operator (18) is one of the form (16) with 
|*(zi, z 2 ,... , z N )) given by 

|#(zi, z 2 , . . . , z N )) = a ZuZ2: ... :ZN I zx, z 2 ,... , z N ) 

+ P Zl ,z 2 ,...,z N **•■■■>** \-Zi, -z 2 ,... , -z N ) , (19) 

where a Zl)Z2) ... jZN and (3 Zl ,z 2 ,... ,z N are real numbers with a 2 Zl Z2 ... ZN + 
P Zl Z2 ... zk = 1- This is so because, as it stands, the state vector (16) 
with 1^(21, z 2 , . . . , Zn)) given by (19) turns out to be the most general 
(pure) state for N spin-1/2 particles and then it will be always possible 
for any eigenvector of the Bell operator (18) to be arranged so as to fit in 
with the form displayed by such Eqs. (16) and (19). Also generally, for 
any eigenvalue A the summation in (16) will extend over all those {nontriv- 
ially) degenerate eigenvectors (19) associated with that given eigenvalue. Of 
course, the coefficients a ZljZ2) ... )ZN and (3 Zl ,z 2 ,... ,z N (or, rather, their prod- 
uct), as well as the phase factors (f> Zl>Z2t ... )ZN will depend on the actual 
value of A. So, for the case in which |A| attains its maximum value 4, (i.e., 
when the inequality is maximally violated), we must have 4> Zl ,z 2 ,...,z N = 

^101 + z 2 (j) 2 H h z N (j) N + nir, and \a zljZ2; ... )ZN \ = \/3 zl>Z2t ... )ZN \ = l/\/2, 

for all z\, z 2 , ■ ■ ■ , Zn (see Eq. (17)). It is easy to show, however, that this 
latter requirement should be fulfilled not only by the maximum possible 
eigenvalue of the relevant Bell operator but also by its largest one. For this 
purpose, let us consider the expectation value of the operator (18) for the 
state vector (16) with |^(-zi, z 2 ,... , zn)} given by (19). This expectation 
value is 

2<X Z 1, Z 2,...,zn (3zi,Z2,...,z N I C ZljZ2i ... jZN \ 

z\,Zi,— , ZN 

x [A a (zi, z 2 , . . . , z N ) + Bp (zi, z 2 ,... , z N ) 
+ C 7 (zi, z 2 ,... , z N ) + D s {zi, z 2 , . . . , z N ) \ , (20) 
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where, for example, C 7 = 7COs(<^ zl)22) ... )2;jv — z\fa[ — Z2§\ — • • • — zjv^jv)- 
Clearly, as the product a Zl ,z 2 ,... ,z N Zl , Z2 ,...,z N factorises out in each of the 
terms involved in Eq. (20), it follows that the largest eigenvalue for i?n must 
fulfil I 

Oi Zl ,z2,— ,z N \ — \flzi,z2,...,z N \ — l/y/2, for all z\, zi,... , zjv- In view 
of Eq. (20) this conclusion holds irrespective of the number of correlation 
functions (A a , Bp, C 7 , Dg, ...) involved. Furthermore, the structure of Eqs. 
(19) and (20) remain unchanged for the general case in which the spin mea- 
surements are carried out along arbitrary directions (in fact, for this case, we 
have only to redefine the meaning of the |zi)'s in (19) as denoting states of 
spin-up (z,i = +1) or -down {z{ = —1) for particle i along some appropriate 
z-axis which, in general, will differ from one particle to the other). Thus, we 
have demonstrated that a state of N spin- 1/2 particles (N > 3) giving the 
largest violation of a Bell inequality must generally consist of a superposi- 
tion of maximally entangled states. This conclusion applies in particular to 
those states giving the maximum possible violation. This is achieved when 
each of the correlation functions attains an appropriate extremum value ±1. 
So, as a direct ("all or nothing") contradiction with local realism arises just 
at the level of perfect correlations, it follows that any state leading to such a 
contradiction should in general involve a superposition of maximally entan- 
gled states. In fact, as we have seen, any state vector yielding the value ±1 
for the correlation function P(fa, fa, ■ ■ ■ , 4>n', must necessarily assume 
the form of Eqs. (16) and (17). 

We conclude by noting that, as expected, this function factorises for a 
general product state of the form |^) = <g) j*!^) <8> ■ ■ ■ <8> \^n), with 
= a Zi \zi) + a- Zi \ - Zi), a Zi = \a Zi \e l6z ^, and \a Zi \ 2 + |a_ z J 2 = 1, 
i = 1, 2, . . . , N. Indeed, by making use of the identity 

2 cos (zi7i + z 2 72 H 1- -2at7at) = 2 Jv cos7i COS72 . . . cos7Ar , 

Z1,Z2,-,Z N (21) 

one can readily see that expression (15) takes the form 

P(fa, fa,... , 4>n; *) = 2^00 COS(01 + r?l) COs{fa +7/2)... C0S{4> N + 7] N ) , 

(22) 

where 77, = 6 Zi - Q- Zi , and clq is a constant with value clq — | a, Zl 1 1 a, Z2 \ ■ ■ ■ | a, ZN | 
x|a_ Zl ||a_2 2 | • • • |a_ Zjv |. Note that ao < 2~ N and then, as it should be, 
\P\ < 1. Obviously, ao reaches its maximum value whenever \a Zi \ = \a- Zi \ = 
l/y/2, for all i. In any case, it is apparent from (22) that for a product state 
the outcome of a spin measurement for any one of the particles becomes 
completely uncorrelated with respect to the outcomes corresponding to the 
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other particles, and then such a state will be unable to yield a violation of 
Bell's inequality. 
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